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Abstract. This paper proposes an efficient tliree fold stratified SIFT matching 
for iris recognition. The objective is to filter wrongly paired conventional SIFT 
matches. In Strata I, the keypoints from gallery and probe iris images are paired 
using traditional SIFT approach. Due to high image similarity at different regions 
of iris there may be some impairments. These are detected and filtered by find- 
ing gradient of paired keypoints in Strata II. Further, the scaling factor of paired 
keypoints is used to remove impairments in Strata III. The pairs retained after 
Strata III are likely to be potential matches for iris recognition. The proposed 
system performs with an accuracy of 96.08% and 97.15% on publicly available 
CASIAV3 and BATH databases respectively. This marks significant improvement 
of accuracy and FAR over the existing SIFT matching for iris. 
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1 Introduction 

Iris is the sphincter having unique flowery random pattern around the pupil. It is an in- 
ternal organ with complex unique features that are stable throughout the lifetime of an 
individual. There has been significant research done in the area of iris recognition using 
global features LL2_J4|. However, these approaches fail to possess invariance to affine 
transformations, occlusion and robustness to unconstrained iris images. Thus, there is 
a stringent requirement to develop iris recognition system suitable for non-cooperative 
images. Keypoint descriptors are invariant to affine transformation as well as partial oc- 
clusion. Scale Invariant Feature Transform (SIFT) is a well known keypoint descriptor 
for object recognition (5). Due to inherent advantages, SIFT is capable of perform- 
ing recognition using non-cooperative iris images [6 1. In SIFT matching approach, the 
difference of Gaussian images are used to identify keypoints at varying scale and ori- 
entation. The orientation is assigned to each detected keypoint and a window is taken 
relative to direction of orientation to find the descriptor vector During recognition, the 
keypoints are detected from gallery and probe images and matching is performed using 
nearest neighbour approach. The challenge with conventional SIFT matching when ap- 
plied to iris recognition is to find texture similarity between same regions of two iris. 
SIFT fails as it does not consider spatial information of the keypoints. To make SIFT 
functional for iris, a novel matching approach has been developed that combines spatial 
information along with local descriptor of each keypoint. 




(a) Strata I: ;/ = 98 (b) Strata II: ;/ = 65 




(c) Strata III: 77 = 54 



Fig. 1. Matches (77) at different strata obtained between two instances of same individual taken 
from CASIAV3 

The organization of the paper is as follows. Conventional SIFT matching approach 
for local descriptors is discussed in Section |2] Proposed stratified SIFT matching is ex- 
plained in detail in Section[3] The job of filtering is done in two sequencial steps: Gra- 
dient based filtering and Scale based filtering. Section |4] provides the results obtained 
using proposed approach. Finally, conclusions are given at the end of the paper. 

2 Conventional SIFT Matching 

The matching algorithm plays a significant role in any biometric system as it acts Uke 
a one way gateway through which only genuine matches (if two images are from same 
subject) will pass and imposter matches (if two images are from different subjects) 
are blocked. In local feature matching, the total number of paired keypoints is used 
to find the authenticity of an individual. Let / be the set of all images available in 
the iris database. For understanding, /,„ be a gallery iris image and /„ be a probe iris 
image where /,„,/„e /. The matching algorithm validates /„ against /,„. In the conven- 
tional SIFT matching, for each keypoint in /,„ the Euclidean distance is found with 
every keypoint in /„. The nearest neighbour approach pairs the keypoint in /,„ with 
keypoint in /„, iff the descriptor distance between the two (after multiplying with a 
threshold) is minimum |5 1. The two keypoints are paired and removed from the set of 
keypoints detected from /„, and /„. This process is iterated for the remaining keypoints 
until any two keypoints can be matched. This approach performs moderately well for 
unconstrained iris recognition |6 1. However, as SIFT determines image similarity using 
128-dimensional local features only, hence it may wrongly pair (impairment) some key- 
points for iris. Thus, the existing approach is modified using two strata which removes 
impaired matches using spatial information of the matching keypoints and contributes 
in achieveing better recognition accuracy. 



3 Stratified SIFT Matching 



In the proposed paper an effort is made to improvise the conventional SIFT matching. 
The pupil and iris circles are assumed to be concentric, hence all localized images have 
pixel size 2r x 2r, where r is the radius of iris. The pupil center as well as iris centre 
are located at (r, r). Therefore the localized images do not have transformation due 
to translation. However, there is a possibility of iris images being transformed due to 
rotation (tilt of subject's head), scaling (change in camera to eye distance) or both ||6l. 
The SIFT matching algorithm matches keypoints that have similarity between the local 
descriptors (as discussed in Section |2]i but fails to conform to spatial relationship. The 
removal (filtering) of impairments by the proposed approach retains only those matches 
that are more probable to be potential matches. Let be the set of m\ keypoints found 
in /„, and K„ be the set of n\ keypoints found in /„ by applying SIFT detector. These 
sets of keypoints are used to comprehend the stratified SIFT matching. 

3.1 Strata I: SIFT Matching 

Let R be the ordered set containing the matches between K„i and K„ by conventional 
SIFT matching as discussed in Section|2] Hence, R contains only those pairs (/, j) where 
keypoint in K,„ is matched with keypoint in K„ as shown in Fig. |l(a)| Let 77 be the 
number of matches found where 77 e [0, min(mi, ni)]. As set R is generated solely on the 
basis of local descriptor property it may wrongly pair keypoints from different regions 
of iris. Hence, there is a need to combine spatial information with local descriptor to 
filter out impaired keypoints as discussed in subsequent strata. 

3.2 Strata II: Gradient Based Filtering 

In this strata gradient based filtering is performed to remove impairments from R. To 
compute gradient for each pair of keypoints (/, j) in R, the angles are obtained from 
respective image centers (r, r). Thus, fl,- is computed from /,„ and (pj is computed from 
/„. The angle of rotation fork''' pair is calculated as - {(pj-©i) mod 360° (depicted in 
Fig. |3(a)| l. Considering SIFT to be completely _^aw/ess (due to robustness property, no 
false match is found) and efficient (due to property of repeatability, all possible matches 
are found) [5J; the value of yu derived should be same VA:. But in practice, SIFT does not 
give such precise matches. Thus, it is difficult to obtain unique value of y even when 
and /„ belong to the same subject. Rather a distribution of y is obtained. A histogram is 
plotted with x-axis comprising bins with a range of values of y, and y-axis comprising 
number of matches falling in a particular bin as shown in Fig. [2] The number of bins in 
the histogram {nobins) is subject to implementation issue. In proposed system, nobins 
is taken as 10. The distribution of y gives a single peak in case the two iris images {!„ 
and /„) are from the same subject. In contrast, no distinct peak should be found in case 
the two iris images are from different subjects. There may be error due to discretization 
of bins, so two adjacent bins of the peak are combined to improve peak density (number 
of matches). The idea is to find whether the density of the peak exceeds the boundary 
criteria. It is inferred that a peak is strong if the density exceeds certain higher bound 
{hp% of total number of matches). 
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Fig. 2. Distribution of y for number of matches between two instances of same individual from 
CASIAV3 



Likewise peak is weak if the density is less than lower bound (lp% of total num- 
ber of matches). If a strong peak is found, an angular range is specified around the 
peak. Those matches in R for which y are not within the angular range are found to be 
impaired and removed from R to generate Rinter. 




(a) y = (0, - e,) mod 360° (b) = | 

Fig. 3. (a) Gradient computation in Strata II, (b) Local scaling factor computation in Strata III 



For example, as shown in histogram in Fig. |2] the peak is found at 0''' bin which 
represents gradient value of 0° to 36° with a central value of 18°. Hence only those 
pairs having angular range between (18 + 90) mod 360° are retained. Thus, it is evident 
that Rinter c R after removing some impairments. Fig jl(b)| shows paired keypoints in 
Rinter with considerable reduction of rf. If no strong peak is found it is inferred that all 
matches in R are faulty, and removed. As a result Rinter becomes empty. 



3.3 Strata III: Scale Based Filtering 



In this strata further fihering of Rinter is performed on the basis of global and local 
scaling factor between the gallery and probe images. The global scaling factor {sf) 
between two images is defined as ratio of probe iris radius (r„) to gallery iris radius (rm). 
A scale range with certain tolerance around sf is empirically taken to handle aliasing 
artifact. From implementation perspective, the scale range is taken as {sf + Q.2). In ideal 
case, if gallery and probe belong to same individual the scaling factor between all paired 
keypoints should be unique. However, this does not hold good in practical scenarios. 

During filtering, for each element in Rinter, two Euclidean distances are calculated- 
(a) d\: distance of keypoint of /,„ from its center and (b) dy. distance of keypoint 
of /„ from its center. Local scaling factor ((/r) for each element of Rinter is calculated 
as (/r = d2ld\ (shown in Fig. |3(b)| i. Matches having if/ within scale range discussed 
above qualifies to be potential and stored in Rnew, else are labeled as faulty and filtered. 
Fig ] l(c)| shows paired keypoints in Rnew after further reduction of 77. 

4 EXPERIMENTAL RESULTS 

The proposed stratified SIFT matching is tested on publicly available BATH |8| and 
CASIAV3 (CASIAV3) |91 databases. Database available from BATH University in- 
cludes images from 50 subjects (20 images per subject from both the eyes). CASIAV3 
(CASIAV3) comprises 249 subjects with total of 2655 images from both the eyes. The 
experiments are carried out on 2.13GHz Intel(R) CPU using Matlab. To validate the 
system performance some standard error measures [71 are usecQ. The results are car- 
ried out in three diff'erent strata as given in Table [1] In Strata I, the two iris images are 
matched using conventional SIFT approach. This approach performs with an accuracy 
of 85.8 1 % on CASIAV3 database. Likewise, for BATH database an accuracy of 97.04% 



Table 1. Performance measures for stratified SIFT matching 



Databases — » 


CASIAV3 


BATH 


Approach J, 


FAR 


FRR 


ACC 


d' 


FAR 


FRR 


ACC 


d' 


Strata I (Conventional SIFT) 


17.48 


10.91 


85.81 


1.20 


1.57 


4.35 


97.04 


2.73 


Strata II (Gradient based Filtering) 


2.49 


5.45 


96.03 


2.46 


0.97 


6.09 


96.47 


2.81 


Strata III (Scale based Filtering) 


2.39 


5.45 


96.08 


2.20 


1.34 


4.35 


97.15 


2.90 



is obtained. To improve the performance of the system, the objective of the proposed 
research is to reduce false acceptances. In Strata II, the impairments are removed us- 
ing gradient filtering which significantly increases the seperability measure between 
false and genuine matches as indicated by d' values given in the Table [1] Further im- 
provement in separability and accuracy are brought by scale filtering in Strata III. The 
accuracy values are plotted against change in number of matches as shown in Fig. |4(a)| 



' FAR: False Acceptance Rate, FRR; False Rejection Rate, ACC: Accuracy, d': d-prime value 




Fig. 4. (a) Accuracy curve for three strata, (b) ROC curve for three strata on CASIAV3 



The Receiver Operating Characteristic (ROC) curves [7 1 for three different strata 
are shown in Fig. |4(b)| The distribution of genuine and imposter scores after Strata 
III is shown in Fig. |5] All graphical results are obtained for CASIAV3 and similar 
observations are made for BATH database. 
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Fig. 5. Final histogram of scores on CASIAV3 



5 CONCLUSIONS 



In this paper, a novel stratified SIFT matching technique is proposed that improvises 
conventional SIFT by removing wrong pairs. This approach provides boost in accu- 
racy due to considerable reduction in FAR. The FAR is reduced by 15.09% and 0.23% 



for CASIAV3 and BATH respectively. From the results it has been observed that the 
proposed algorithm is completely flawless, i.e., matches removed are guaranteed to be 
wrong matches whereas it is not completely efficient, i.e., all impairments by SIFT are 
not guaranteed to be filtered. However, the gain in accuracy is substantial which marks 
its applicability for unconstrained iris recognition. 
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